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Final  Technical  Report 


Computational  Nonlinear  Optics: 
Femtosecond  Atmospheric  Light  String  Applications 

Grant  #:  AFOSR  F49620-03-1-0194 

Grant  Period  :  January  1,  2003  -  December  31,  2005.  No-cost 
extension  to  June  29,  2006. 

Principal  Investigator  :  J.V.  Moloney 

Project  Goals 

The  key  goal  of  this  project  was  to  provide  a  rigorous  theoretical  and  computational 
framework  to  check  the  feasibility  of  various  applications  in  atmospheric  light  string  and 
other  nonlinear  optical  propagation  phenomena  requiring  resolution  of  extremely  fine 
features.  Fully  “optical  carrier  resolved”  vector  femtosecond  pulse  propagators  will 
provide  the  basis  for  investigating  a  broad  landscape  of  nonlinear  optical  interactions 
including:  white- light  super- continuum  generation  due  to  electromagnetic  optical  shocks, 
directionally  narrow  white- light  scattering  from  water  droplets,  backward  Raman 
generation.  Research  would  be  carried  out  into  deriving  more  realistic  physical  material 
interaction  models  appropriate  to  femtosecond -duration  pulse  interactions. 

The  project  proposed  a  multi-disciplinary  collaboration  with  the  experimental  group  of 
Professor  J.C.  Diels  at  the  University  of  New  Mexico,  with  Dr.  James  Murray  of  Lite 
Cycles,  Tucson  and  with  Professor  FHM  Faisal  and  Dr.  Andreas  Becker  at  the  University 
of  Bielefeld,  Germany.  The  collaboration  with  the  experimental  group  of  Professor  Diels 
was  to  provide  an  experimental  validation  route  to  the  modeling.  Dr.  Murray  acted  as  a 
project  consultant  providing  input  on  potential  Air  Force  applications.  Professor  Faisal 
and  Dr.  Becker  are  world  experts  on  ab  initio  calculations  of  ultra-  intense  laser  pulse 
interactions  with  molecules,  including  N2,  and  O2,  key  constituents  of  air. 

Supplemental  funding  was  added  to  the  original  grant  to  carry  out  a  one -year  study  of  the 
feasibility  of  selectively  ablating  aircraft  engine  turbines  as  a  means  of  healing 
microscopic  cracks.  A  technical  report  “Ultrafast  Laser  Pulse  Machining  of 
Coated  Nickel- Base  Alloy  “  was  submitted  to  the  Air  Force  on  completion  of  this  project 
in  2004. 

A  parallel  effort  was  initialized  in  computational  nanophotonics  in  anticipation  of  the 
need  to  eventually  interface  the  pulse  propagators  to  Maxwell  FDTD  solvers.  Uniform/ 
non-uniform  mesh  and  an  adaptive  space  and  time  3D  FDTD  solvers  were  developed.  A 


goal  here  was  to  better  understand  ultrashort  pulse  interaction  with  dielectric  disks  and 
spheres  as  a  prelude  to  studying  water  droplets. 

Executive  Summary 

The  continued  success  enabled  by  this  grant  maintains  researchers  at  the  Arizona  Center 
for  Mathematical  Sciences  (ACMS)  at  the  forefront  worldwide  in  modeling  ultra-short 
intense  laser  pulse  propagation  in  various  linear  and  nonlinear  media.  The  large  number 
of  publications  in  Physical  Review  Letters  and  other  top  quality  journals  as  a 
consequence  of  the  research  generated  under  is  project,  is  testament  to  the  high  regard  of 
the  work  in  the  eyes  of  the  international  community.  The  key  breakthrough  result  was 
the  publication  of  the  paper  “Unidirectional  pulse  propagation  equation”  by  M.  Kolesik, 
J.V.  Moloney  and  M.  Mlejnek  in  Physical  Review  Letters,  89,  2839021  (2002).  This 
work  on  UPPE  laid  the  foundation  for  the  many  successful  achievements  of  the  present 
grant  and  provided  the  first,  fully-carrier  resolved,  unidirectional  vector  pulse  propagator 
that  could  evolve  3D  arbitrarily  shaped  laser  pulses  over  many  meter  distances. 

Moreover,  this  UPPE  model  provided  the  first  systematic  means  of  transitioning 
seamlessly  to  the  famous  Nonlinear  Schrodinger  Equation  pulse  propagation  model  and 
to  the  many  more  recent  propagation  models  proposed  as  an  improvement  on  NLSE.  In 
stark  contrast,  3D  FDTD  Maxwell  solvers  are  so  computationally  intensive  that  they  are 
restricted  to  fine  three-dimensional  domains  spanning  typical  computational  volumes  of 
dimension  lOOpmxlOOpmx  1000pm.  FDTD  vector  Maxwell  solvers  are  further  restricted 
in  application  by  the  rather  simplistic,  spectrally  local,  single  and  multiple  pole  Lorenz 
light- matter  models.  United  States  Air  Force  Captain  Thomas  Niday,  completed  his  Ph.D 
studies  in  Optical  Sciences  at  the  University  of  Arizona  working  directly  on  this  project. 
The  title  of  his  Ph.D  Thesis  is  entitled  “Stability  and  Transient  Effects  in  Ultrashort  UV 
Filaments”.  Dr.  Niday  successfully  defended  his  Ph.D  on  July  15,  2005. 

The  first  significant  application  of  UPPE  under  this  project  was  to  atmospheric  light 
string  propagation.  The  complex  nonlinear  physics  occurring  in  the  intense  focal  region 
of  the  critically- self- focused  pulse  (white- light  supercontinuum  generation  spanning  the 
entire  UV  -  IR  spectrum,  optical  breakdown  of  air  with  plasma  generation  that  blue-shifts 
the  supercontinuum  and  emits  an  intense  burst  of  THz  short  range  radiation  from  the 
recombining  plasma)  necessitates  a  nonlinear  pulse  propagator  that  can  resolve  light- 
matter  coupling  over  these  huge  spectral  ranges.  The  experimental  team  of  Professor 
Diels  at  UNM  were  tasked  with  setting  up  light  string  experiments  in  the  UV.  Because  of 
technical  details  and  lack  of  specialized  hardware  resources,  this  effort  was  not 
productive  and  terminated  after  one  year.  Dr.  James  Murray  remained  as  consultant  on 
the  project  in  spite  of  his  moving  from  Lite  Cycles  to  Arete  Associates  in  2004. 
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In  this  AFOSR  funded  project  “ Computational  Nonlinear  Optics:  Femto  Atmospheric 
Light  String  Applications"  AFOSR  F49620-03-1-0194,  we  further  developed  and 
demonstrated  a  novel  ultrashort  unidirectional,  optical  carrier-resolved,  vector  pulse 
propagator  (UPPE).  The  development  of  this  model  was  necessary  because  of  the 
obvious  shortcomings  of  classical  nonlinear  envelope  pulse  propagation  models  such  as 
Nonlinear  Schrodinger  Equation  (NLSE)  and  its  many  variants  when  applied  to 
ultrashort,  high  intensity  pulse  propagation  phenomena  in  bulk  media.  The  NLSE  is  a 
slowly  varying  envelope  propagator  for  near  paraxial  pulses  that  necessarily  is  restricted 
to  spectrally  narrow  nonlinear  excitations.  A  detailed  derivation  of  NLSE  with  higher 
order  corrections,  as  an  asymptotic  expansion  in  a  small  parameter  can  be  found  in  the 
textbook  Nonlinear  Optics  by  Moloney  and  Newell  [Moloney 2004],  It  was  recognized 
that  newly  emerging  high-  field  ultrashort  pulse  sources  involve  both  strong  deviation 
from  paraxiality  and  ultrabroadband  excitation.  Consequently,  systematically  adding 
additional  terms  to  NLSE  in  a  self-consistent  fashion  should  eventually  fail. 

The  UPPE  can  be  formulated  in  both  a  time-  and  z-propagated  form.  For  the  present 
purposes  we  present  just  the  full  vector  z-propagated  version  that  connects  most  naturally 
to  the  many  nonlinear  envelope  models  presented  in  the  recent  literature. 

(<A  z)  =  ikzA\  (<*>»  z)+  Z  * 

s=\,2 

where 

K  =yl(02£(co)/c2 -k]-k) 

is  the  plane -wave  propagation  constant  for  angular  frequency  ro  and  transverse 
wavenumbers  kx,  ky.  This  z-propagated  form  of  UPPE  describes  the  evolution  of  spatio- 

temporal  spectral  amplitudes  Ek  k  (to,  z)  of  the  optical  field  along  the  propagation 

coordinate.  The  real  total  field  amplitude  E(x  ,y,z,  t) )  is  obtained  by  Fourier 

transforming  the  spectral  amplitude  at  a  fixed  space  point  z.  This  UPPE  has  all  of  the 
linear  absorption/dispersion  properties  lumped  into  the  first  term  on  the  right  and 
nonlinear  polarization  and  current  sources  constitute  the  remaining  terms.  This  equation 
encodes  the  complete  linear  optical  response  (absorption/dispersion)  in  the  leading  term 
on  the  right  and  can  handle  arbitrarily  strong  departure  from  paraxial  behavior  (extreme 
linear  or  nonlinear  focusing  down  to  ~X,  for  example).  The  nonlinear  source  terms  are 
not  discussed  in  detail  here  as  it  is  the  physically  self-consistent  formulation  of  these  that 
represents  a  major  thrust  of  the  proposed  effort. 
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For  many  practical  applications  the  scalar  version  of  this  equation  suffices 
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In  this  form,  we  can  seamlessly  derive  the  many  nonlinear  envelope  models  of  the  recent 
literature  [Kolesik2004UPPE]. 


As  an  illustration  we  will  derive  NLSE  and  a  recent  nonlinear  envelope  equation  (NEE) 
proposed  by  Brabec  and  Krausz  in  1997[Brabecl997].  We  assume  no  free  charges  and 
drop  the  last  current  source  term.  Writing  z-UPPE  in  a  compact  form 


dzEkxx  ("’z)  =  iKEkk  z)+iQpk„k,  (">z) 


where  K  (kx  Jcy  ,<wj  =  ^<w2e(  aty/c2  -  k^  -  l?y  is  the  linear  field  propagator  in  the  spectral 
representation  and 
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will  be  called  the  nonlinear  coupling  term.  We  will  have  to  make  some  simplifying 
assumption  about  the  polarization  in  what  follows  in  order  to  connect  to  familiar 
envelope  models. 


To  obtain  envelope  equations  we  must  express  the  electric  and  polarization  fields  in 
terms  of  envelopes  by  factoring  out  the  carrier  wave  at  a  chosen  reference  angular 
frequency  (0R  with  corresponding  wavenumber  kR  =  K(0,0,O)r)  : 

E(x,y,z,t)  =  A(x,y,  z,  t)  d(kR2-mR,) 
and  similarly  for  the  polarization  P. 

The  procedure  will  be  to  Taylor  expand  the  coefficients  K  and  Q  above  and  use  the 
standard  identities: 


(to-  (0R )  o  id,  ikx  <=>  dx  iky  <=>  dy  3.0  ik(coR)  +  d.  . 

The  classical  NLSE  equation  follows  by  expanding  the  linear  term  K  to  second  order  in 
(ty-CDfl)  and  kx,ky: 

K ( kx ,ky , (o)  =  Jco2  e(  aj) / c2  -  kj  -  k2y 

and  ignoring  completely  the  frequency  dependence  of  the  nonlinear  coupling  term  Q  i.e 
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In  order  that  we  recover  the  classical  NLSE  we  make  a  simple  assumption  for  the 
nonlinear  polarization  i.e  assume  an  optical  Kerr  nonlinearity 

P  =  2e(/i(coR)  n2l A 

Inserting  these  approximations  into  z-UPPE  we  obtain 


The  physical  restrictions  of  NLSE  are  evident  from  this  approximation  procedure. 
Approximating  K  to  second  order  in  frequency  and  transverse  wavenumber  amounts  to 
the  paraxial  and  quasi- monochromatic  approximation  for  linear  wave  propagation.  The 
approximation  that  Q  is  a  constant  also  requires  a  spectrally  local  interaction. 

The  nonlinear  envelope  model  (NEE)  derived  by  Brabec  and  Krausz  in  1997  is  widely 
used  in  the  modeling  of  high  power  ultrashort  pulse  propagation  problems.  For 
completeness  and  as  a  contrast  to  NLSE,  we  derive  this  NEE  model  now. 

We  again  approximate  the  linear  propagator  K  by  its  paraxial  version: 


This  is  essentially  the  same  second  order  (paraxial)  Taylor  expansion  in  transverse 
wavenumber  with  some  minor  approximations.  Namely,  we  replace  n(coi)  -+n(coR)  in 

the  denominator  of  the  diffraction  term  and  thus,  partly  neglect  chromatic  dispersion.  We 
retain  the  full  linear  dispersion  and  for  convenience  rewrite  K  so  as  to  explicitly  extract 
the  group  velocity  term 
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contains  all  higher  order  temporal  dispersion  terms. 

We  next  approximate  the  nonlinear  coupling  term  Q  as: 
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Putting  these  approximations  back  into  z-UPPE  we  obtain 
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Finally,  transforming  to  the  real  space  representation 
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This  is  the  NEE  equation  used  extensively  in  ultrashort  pulse  propagation  problems.  It 
has  the  obvious  advantage  of  including  temporal  dispersion  to  all  orders  but  still  ignores 
non-paraxial  effects.  Notice  the  appearance  of  the  inverse  operator  coupling  time  and 
space  -  this  is  often  truncated  using  a  binomial  expansion  with  potential  disastrous 
consequences  for  the  study  of  supercontinuum  generation. 

The  class  of  problems  that  we  intend  to  tackle  in  this  proposal  will  generally  involve 
strong  non-paraxial,  ultra- broadband  propagation  phenomena  and  truncated  models  such 
as  the  above  are  not  relevant.  Moreover,  nothing  has  been  said  about  how  one 
approximates  the  induced  nonlinear  polarization  or  current  source  terms.  Approximations 
made  in  the  literature  for  these  terms  are  even  more  severe  and  deriving  better  nonlinear 
source  models  represents  a  major  proposal  thrust. 


II.  Success  Stories  of  UPPE 

We  now  briefly  encapsulate  some  of  the  successes  achieved  to  date  using  the  UPPE 
propagator  model.  This  model  has  been  applied  to  femtosecond  atmospheric  light  string 
propagation,  simultaneous  generation  of  supercontinuum  and  third  harmonics  in  air  and 
water,  generation  of  nonlinear  X-waves  in  normally  dispersive  and  O- waves  in 
anomalously  dispersive  condensed  media  and  a  reduced  effective  1 D  model  to  study 


supercontinuum  shaping  in  sub-micron  and  photonic  crystal  fiber  cores.  These  successes 
with  rather  naive  nonlinear  dispersion  models  can  be  ascribed  to  the  fortuitous  dominance 
of  the  linear  dispersion  in  each  scenario.  A  goal  of  the  present  proposal  will  be  to  derive 
computationally  feasible  and  physically  self-consistent  models  that  capture  the  full  linear 
and  nonlinear  dispersive  and  absorptive  properties  of  the  relevant  material  system. 


11.1  Femtosecond  Atmospheric  Light  Strings 


Our  group  at  the  Arizona  Center  for  Mathematical  Sciences  (ACMS)  has  played  the  lead 
internationally  in  establishing  a  robust  theory  and  modeling  of  femtosecond  light  string 
propagation  in  air.  The  mechanism  of  nonlinear  spatial  replenishment  was  proposed  by  us 
in  1998  [Mlejnekl998]  and  shortly  thereafter,  experimentally  confirmed  by  the  group  in 
Laval,  Quebec.  It  was  proposed  by  us  at  the  time  that  the  persistent  nonlinear  chaotic 
generation  of  intense  light  strings,  accompanied  by  plasma  channels  could  be  identified 
with  the  broad  pulse  background  acting  as  a  reservoir  that  feeds  individual  critical  self- 
focusing  events  within  the  much  broader  pulse.  While  the  initial  pulse  waist  is  on  the 
order  of  centimeters,  individual  light  strings  have  diameters  on  the  order  of  100  pm  and 
the  dilute  plasma  channels  generated  by  the  light  strings  even  smaller  50  -  80  pm 
diameters.  Our  initial  model  for  light  string  generation  was  a  generalized  NLSE  equation 
coupled  to  a  Driide  plasma  model.  This  model  proved  important  in  isolating  the  role  of 
various  physical  mechanisms  operative  during  the  initial  nonlinear  interaction  over  a  few 
meters.  We  were  very  much  aware  at  the  time  that  this  NLSE  model  had  some  serious 
drawbacks.  For  example,  its  narrowband  (slowly  varying  envelope)  limitation  raised 
concerns  about  properly  describing  supercontinuum  generation  and  we  moved  on  to 
deriving  UPPE.  Of  equal  concern  however  was  that  we  were  describing  the  nonlinear 
interaction  terms  by  fixed 
wavelength  pre-computed 
optical  Kerr  coefficients  and 
multiphoton  cross-sections  - 
yet  the  ultra- Broadband 
generated  supercontinuum 
has  spectral  intensities  in  the 
10’s  of  GW  range  that  are 
fully  consistent  with  strong 
nonlinear  coupling.  UPPE 
solved  the  former  problem 
but  the  latter  remains  open  as 
we  need  quantum  mechanical 
models  of  the  light- matter 
interaction  -  the  latter  are  as 
computationally  as  intensive 
or  more  so  than  the  EM 

propagation  equation.  One  fundamental  question  that  has  not  been  resolved  in 
experimental  femtosecond  light  string  propagation  studies  is  the  appearance  of  a 
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Figure  2  Real  water  susceptibility  and  an  artificially 
constructed  susceptibility  that  departs  for  the  real  at  high 
frequencies  well  outside  the  initial  pulse  spectral  bandwidth. 


relatively  quiescent  regime  with  weak  plasma  generation  subsequent  to  the  initial  violent 
critical  collapses.  Experiment  indicates  that  the  light  strings  seem  to  stay  collimated  even 
though  there  appears  to  be  no  indicator  of  plasma  generation  i.e  detectable  THz  signals. 
This  is  consistent  with  some  form  of  saturation  of  the  optical  Kerr  effect.  We  have  very 
preliminary  evidence  from  full  quantum  mechanical  solutions  of  hydrogen-  like  atoms 
that  such  a  saturation  effect  is  real  [Kano2006].  Hence,  although  many  of  the  qualitative 
features  can  be  understood  with  realistic  linear  dispersion  and  non-dispersive  nonlinear 
effects,  this  latter  important  observation  must  be  ascribable  to  a  nonlinear  dispersive 
response. 

Using  UPPE  we  were  able  to  show  that  the  spectral  extension  of  super- continuum 
generation  about  the  reference  pulse  frequency  is  a  very  sensitive  function  of  the  linear 
material  dispersion  at  spectral  locations  far  removed  from  the  initial  pulse  spectrum.  We 
performed  comparative  simulations  in  water  and  gases  to  reveal  the  mechanism  that 
controls  the  supercontinuum  extent  and  its  dependence  on  the  medium  band -gap 
[Kolesik2003PRL:  Kolesik2003APB]. 

Figure  2  shows  two  linear  material  susceptibilities  for  water  and  an  artificial,  water- like 
medium.  These  were  chosen  such  that  the  bulk  of  the  SC-generating  pulse  propagated 
essentially  unchanged.  Figure  3  shows  on  the  left  how  the  spectral  broadening  depends 
on  input  pulse  energy.  On  the  right,  we  show  that  the  modification  of  the  real  part  of  the 
dielectric  susceptibility  far  away  from  the  influence  of  the  exciting  pulse  bandwidth  can 
significantly  alter  the  spreading  of  the  generated  SC.  This  implies  a  strong  sensitivity  to 
the  linear  dispersion  landscape. 

These  simulations  allowed  us  to  correct  and  extend  the  so  far  accepted  scenario  for  SC 
generation  and  its  dependence  on  the  medium  band- gap.  An  effective  three- wave  mixing 
paradigm,  [Kolesik2004PRL]  that  can  be  used  to  explain  qualitatively  a  number  of 
effects  in  femtosecond  pulse  propagation  in  general,  also  emerged  as  a  product  of  our 
qualitative  simulation  studies. 
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Figure  3  Left:  Measured  spectral  superbroadening  in  water  as  a  function  of  increased  pulse  energy.  Right: 
Spectral  extent  of  the  generated  supercontinuum  with  and  without  the  artificial  water  susceptibility  response. 
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11.2  Supercontinuum  and  Third  Harmonic  Generation 


More  recently  we  succeeded  in  reproducing  experimental  results  showing  simultaneous 
SC  and  THG  generation  in  air  under  weak  and  strong  linear  focusing  conditions.  By 
using  a  short  focal  length  lens,  all  nonlinear  effects  occur  near  the  localized  high  intensity 
focal  spot  making  nonlinear  propagation  effects  less  important.  Consequently,  the 
generated  SC  was  localized  in  frequency  and  both  fundamental  and  THG  spectral 
components  remained  essentially  decoupled  during  propagation.  These  results  are  shown 
in  Figure  4. 

The  very  strong  intermixing  of  fundamental  and  third  harmonic  components  on  the  right 
in  Figure  4  makes  it  unlikely  that  a  classical  slowly- varying  individual  FF  and  TH 
envelope  field  approach  will  work. 

Our  simulations  were  the  first  to  include  TH  and  SC  generation  in  a  unified,  self- 
consistent  way  [Kolesik2006].  Our  results  make  it  possible  to  correct  some  previously 
published  conclusions.  In  particular,  we  have  shown  that  the  SC  generation  mechanism  is 
independent  of  the  TH  generation  mechanism,  though  both  effects  are  controlled  by  the 
same  underlying  effective  three- wave  mixing  processes  [Kolesik2004PRL].  We  have  also 
shown  that  it  is  possible  to  clearly  distinguish  experimentally  the  origin  of  the  high- 
frequency  radiation  in  the  far- field  spectrum  even  when  the  TH  and  SC  components 
appear  to  merge  as  illustrated  in  Figure  4. 
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Figure  4  Left:  On-  and  off-axis  emission  spectra  for  a  tightly  focused  intense  femtosecond  pulse  in  water. 
The  fundamental  (FF)  and  third  harmonic  (TH)  don’t  overlap  spectrally.  Right:  Strongly  admixed  FF  and 
Tl  I  spectral  components  due  to  strong  SC  generation  subsequent  to  critical  collapse  in  weakly  focused 
geometry. 


11.3  Nonlinear  X-  and  O-Waves  in  Water 


The  subtle  interplay  between  critical  collapse  (self-  focusing)  and  normal/anomalous 
dispersion  in  condensed  media  has  been  spectacularly  demonstrated  in  a  series  of 
important  experiments  by  a  group  in  Como,  Italy.  By  loosely  focusing  an  intense 
femtosecond  pulse  in  water,  they  observed  a  very  strong  X- feature  in  the  spectrally 
resolved  far- field  of  the  trapped  pulse  when  the  carrier  wavelength  was  such  that  the 
material  exhibited  normal  dispersion.  More  recently,  nonlinear  O- waves  were  observed 
by  the  same  group  when  the  central  pulse  wavelength  was  chosen  to  correspond  to 
anomalous  material  dispersion.  These  propagation  phenomena  are  illustrations  of 
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Figure  5  The  z-UPPE  takes  as  input  the  entire  absorption  and  dis£krsion  spectral  landscape  shown  on  the 
left  and  right,  respectively.  The  numerical  algorithm  senses  all  physical  values  over  the  spectral  range 


strongly  non-paraxial  propagation  requiring  accurate  resolution  of  off-axis  propagating 
components  of  the  propagating  pulse.  The  measured  input  absorption  and  dispersion  data 
for  water  are  shown  in  Figure  5.  Our  work  showed  that  the  propagating  pulse  exhibited 
extremely  complicated  spatio- temporally  chaotic  bursts  giving  the  illusion  of  a  stationary 
self-trapped  nonlinear  waveguide  in  water.  What  is  remarkable  is  that  the  spectrally- 
resolved  far- field  exhibits  a  pronounced  X- feature  implying  the  existence  of  a  universally 


Figure  6  Left:  Spectrally  resolved  far-field  for  a  pulse  where  the  initial  and  generated  SC  is  located  entirely 
within  the  normally  dispersive  region.  Right:  The  pulse  spectrum  lies  in  the  anomalously  dispersive  region  but 
the  generated  spectrum  extends  into  the  normally  dispersive  regime. 


attracting  asymptotic  set  in  the  infinite  degree  of  freedom  problem.  Figure  6  shows  a 
graph  of  the  X-wave  and  O-wave  feature.  The  O-wave  is  polluted  by  some  X-wave 
character  in  the  right  picture  because  the  generated  SC  extends  across  into  the  normally 
dispersive  regime. 

11.4  Self-healing  filaments  and  the  role  of  the  low-intensity 
background 

Computer  simulations  contributed  crucially  to  the  present  understanding  of  femtosecond 
light  strings  in  air  and  other  gaseous  media  as  well  as  in  condensed  matter.  One 
particularly  important  observation,  made  for  the  first  time  by  our  group  [Mlejnekl999] 
several  years  ago,  is  that  the  high-  intensity  light  strings  exchange  their  energy  with  the 
reservoir  of  the  low- intensity  background  of  the  pulse.  Thus,  the  low- intensity 
background  provides  the  mechanism  necessary  for  multiple  hot  filaments  to  survive  over 
large  distances. 

More  recently,  an  experiment  [Courvoisier2003]  suggested  that  even  a  single  filament  is 
extremely  robust  with  respect  to  its  interaction  with  water  droplets  (for  example  in 
clouds).  We  have  applied  our  UPPE  simulator  to  reveal  for  the  first  time  that  it  is  actually 
the  low- intensity  background,  forming  a  pedestal  surrounding  the  central  hot  spot,  that  is 
the  cause  for  the  filament  robustness.  We  have  shown  that  filaments,  after  being 
temporarily  destroyed  in  the  center  by  collision  with  a  water  droplet,  heal  themselves 
over  very  short  distances.  This  self-healing  process  is  driven  by  the  optical  energy  stored 
in  the  low- intensity  background  [Kolesik2004SH].  Note  that  similar  results  were  later 
published  independently  by  another  group  [Skupin2004], 

This  insight,  brought  by  simulation,  provided  the  theoretical  basis  for  understanding  a 
subsequent  open-air  experiment  [Mechain2005]  demonstrating  that  femtosecond  light 
strings  can  propagate  even  in  clouds  and  rain. 
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Figure  7  Left:  Calculated  group  velocity  dispersion  (GVD)  for  silica  fibers  with  two  different  core  diameters: 
a=l.20  um  and  a=0.325  urn.  Right:  Generated  white  lieht  continuum  using  an  NLS  model  and  the  UPPE  model. 


The  role  of  the  low- intensity  background  has  since  been  corroborated  in  several 
experimental  and  theoretical  studies  [Dubietis2004:  Liu2005:  Ackermann2006: 
Gaizauskas2006]. 


11.4  Shaping  the  super-continuum  in  sub-micron  and  photonic  crystal 
fibers 

This  application  of  UPPE  takes  advantage  of  the  very  accurate  representation  of  white- 
light  super-continuum  over  huge  spectral  windows.  The  problem  can  be  reduced  to  a  ID 
simulation  due  to  the  strong  confinement  of  the  light  within  the  fiber  core.  Figure  7 
shows  an  example  of  how  the  fiber  group  velocity  can  be  engineered  by  tapering  a 
conventional  silica  fiber  such  that  the  fiber  core  is  reduced  to  a  sub- micron  diameter.  The 
resulting  GVD  landscape  is  significantly  modified  with  two  zeroes  of  GVD  occurring  for 
the  smaller  core  fiber  [see  Figure  7  Left].  The  right  hand  side  of  this  figure  shows  a 
comparison  of  the  white- light  super-continuum  obtained  by  using  a  generalized  NLSE 
model  (that  includes  accurate  chromatic  and  waveguide  dispersion)  and  the  present  UPPE 
model. 
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